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Abstract 

Phase measurements on a single-mode 
radiation field are examined from a 
system-theoretic viewpoint. Quantum 
estimation theory is used to establish 
the primacy of the Susskind-Glogower 
(SG) phase operator; its phase eigen- 
kets generate the probability operator 
measure (POM) for maximum-likelihood 
phase estimation. A commuting observ- 
ables description for the SG-POM on a 
signal x apparatus state space is derived. 
It is analogous to the signal-band x image- 
band formulation for optical heterodyne 
detection. Because heterodyning realizes 
the annihilation operator POM, this anal- 
ogy may help realize the SG-POM. The 
wave function representation associated 
with the SG-POM is then used to prove 
the duality between the phase measure- 
ment and the number operator measure- 
ment, from which a number-phase uncer- 
tainty principle is obtained, via Fourier 
theory, without recourse to linearization. 
Fourier theory is also employed to estab- 
lish the principle of number-ket causal- 
ity, leading to a Paley-Wiener condi- 
tion that must be satisfied by the phase- 
measurement probability density function 
(PDF) for a single-mode field in an ar- 
bitrary quantum state. Finally, a two- 
mode phase measurement is shown to af- 
ford phase-conjugate quantum communi- 
cation at zero error probability with finite 
average photon number. Application of 
this construct to interferometric precision 
measurements is briefly discussed. 


1. INTRODUCTION 

A classical single-mode radiation field 
is characterized by a spatial-mode pattern 

f (r ), an oscillation frequency u in rad/s, 
and a c- number phasor a. The latter spec- 
ifies both the energy and the initial phase 


shift of the field — we can take H — Nhu) 
to be the mode energy, and <f> to be the 
mode phase, where 

O = y/Ne'*, (1) 

is the polar decomposition of a. When 
the single-mode classical field is quantized, 
its mode pattern and frequency are un- 
changed, but a is replaced by the anni- 
hilation operator a. The phase problem, 
for this single- mode quantum field, has 
long been taken to mean finding a satis- 
factory quantum version of Eq. I. 1 How- 
ever, owing to the noncommutative nature 
of the quantum theory’s operator algebra, 
no such decomposition exists, i.e., there is 

no observable $ such that 2 

a = yftie*. ( 2 ) 

One may quibble about the order of the 
amplitude and phase terms on the right- 
side of Eq. 2, or prefer the use of TV + 
/ = aa^ in lieu of TV = a^d, etc., but the 
essential issue is the nonexistence of the 
observable <f>. 

Until recently the Susskind-Glogower 
(SG) phase operator, 2 

e'* = (dat)~ 1//2 d, (3) 

has seemed to provide the best quantum 
description of phase. The SG operator is 
non-Hermitian, and its quadratures, 

cos(<£) = Re(e i< ^ > ), and (4) 

s\n(4>) = Im(e^), (5) 

are noncommuting observables which fail 
certain reasonable conditions that the co- 
sine and sine of a phase should meet. For 
example, it turns out that 

(i/>|cos(<£) 2 |t/>) + <^|sin(0) | rp) < 1, (6) 
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unless the state \tp) is orthogonal to the 
vacuum, |0). On the other hand, the SG- 
based commutator 

[/V, sin (<£)] = icos(<£), (7) 

does lead to the oft-employed number 
phase uncertainty principle, ANA<f> > 
1/2, under a high-mean-field linearization. 

Lately, there has been intense renewed 
interest in the quantum phase problem. In 
what follows we will review some of the 
recent quantum phase work of Shapiro, 
Shepard, and Wong, 3, 4 and present some 
new results. Because the effort of Shapiro 
et al. originates from a quantum esti- 
mation theory viewpoint, that tack will 
taken here as well. Because the formal- 
ism of Shapiro et al. relies on the prob- 
ability operator measure (POM) descrip- 
tion of quantum measurement — a general- 
ization of observables not well known in 
the physics literature — we will begin with 
a brief tutorial on POM’s. 


2. POM REVIEW 

The textbook approach to quantum 

measurement is through observables . 5 For 
example, consider the quadrature compo- 
nents of the single- mode field’s annihila- 
tion operator, i.e., 

di = Re(d), and (8) 

a 2 = Im(a). (9) 

These are continuous-spectrum observ- 

ables. In other words, they are Hermitian 
operators 


a] = aj, for j = 1,2, (10) 


Measurement of a quadrature operator, 
when the system is in state gives a 
continuous-valued, classical random vari- 
able with PDF 

P( a j 11^)) = lj(«j|'0)| 2 , 

for -oo < atj < oo, j = 1,2. (14) 

For this classical probability density to be 
correct, for all possible |^), it must satisfy 

p( I IV>) ) > o, 

for -oo < Qtj < oo, ji — 1,2, (15) 


and 


f°° 

/ da jP ( aj IIV')) = l, 

for j = 1,2. (16) 

These conditions are ensured by Eq. 13, 
which leads to the familiar quadrature 
representations — essentially the position 
and momentum wave functions — given by 


ivo = 


and 

ivo 


/ oo 

da i |ai)n(ai|V») 

-OO 

f 00 

/ dai V’(ai)lai)!, (17) 

J — OO 

/ OO 

da 2 |a 2 )22(<*2|^) 

-OO 

/ da 2 $(a 2 )|a 2 )2, (18) 

J — OO 


with the obvious identifications for ip(ai) 
and ^(a 2 ). Of course, the quadratures are 
noncommuting observables, 


(£ 1 . 02 ] 



(19) 


with complete orthonormal (CON) eigen- 
kets, 


so they cannot be measured simultane- 
ously . 5 


I a j)j — a j I a j )j > 

for —00 < aj < 00,(11) 

= Siatf-aj), ( 12 ) 

I = f datj \aj)jj {a-j | ,(13) 

J — OO 

where I is the identity operator, and £(•) 
is the Dirac delta function. 


The preceding review demonstrates that 
the full specification of observables, i.e., 
Hermitian operators with CON eigenkets, 
is not needed to produce a consistent 
statistical characterization of a quantum 
measurement. For an arbitrary quan- 
tum state, a resolution of the identity — an 
outer-product sum like Eq. 13 — generates 
a proper classical-probability description 
of a quantum measurement. This is the 
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essence of the POM concept. Our princi- 
pal purpose for introducing POM’s is to 
accommodate measurements that are not 
observables on the state space, W, of the 
d-mode. The best way to introduce such 
nonobservable POM’s is through an exam- 
ple. It is well known that the annihila- 
tion operator, a, is not an observable — it 
is non-Hermitian, 

a* / a. (20) 

Furthermore, its real and imaginary 
parts, di and a 2 , are noncommuting 
observables — a cannot be measured in the 
usual textbook sense. However, the anni- 
hilation operator does have eigenkets — the 
coherent states, 6 

d\a) = a\a) } 

for a € C, (21) 

where C is the complex plane. These states 
are not orthonormal, i.e., 

(<*» - 

exp(-^|a'| 2 - ^|a| 2 +a'*a) ,(22) 

which is a consequence of the nonvanishing 
commutator 

[a, a+] = /, (23) 

implied by Eq. 19. Nevertheless, the 
coherent states are complete, in fact 
N ouercomplete — they form a resolution of 
the identity 

I = -[ d 2 a|a)R (24) 
7T Ja€C 


hold, for all |^)- 

The preceding POM has long been 
known. 7, 8| 9 It represents a measure- 
ment, in the POM sense, of the annihi- 
lation operator a: the d-eigenkets gener- 
ate the measurement statistics, and the a- 
eigenvalues are the resulting observation 
values. This parallels the usual observ- 
ables description: the observable’s eigen- 
kets generate the measurement statistics, 
and the associated eigenvalues are the ob- 
servation values, cf. Eq. 14. 

The POM formulation is not in conflict 
with the conventional dictum that only ob- 
servables can be measured. Any nonob- 
servable POM on H can be represented as 
a collection of commuting observables on 
some larger state space which describes 
the original system interacting with an 
appropriate apparatus. 7 The most famil- 
iar example of this genre is optical het- 
erodyne detection of a single-mode signal 
field, which provides both a commuting- 
observables description and a physical re- 
alization for the d-POM. 

In optical heterodyne detection, 8, 9 > 10 
a signal field of frequency v is mixed with 
a strong local-bscillator (LO) field of fre- 
quency v - vi p on the surface of a pho- 
todetector. With a unity-quantum ef- 
ficiency detector, and two-channel lock- 
in amplification at the intermediate fre- 
quency (i>7f)j this arrangement produces 
a complex-valued, classical random vari- 
able, j/, whose measurement statistics are 
identical to those of the operator 9 

y = (ds <8> • (28) 


which defines the d-POM. The outcome of 
the cLPOM is a complex- valued, continu- 
ous classical random variable with PDF 

p(a\\rl>)) = -|(q|^)| 2 , 

7T 

for a € C, (25) 

when the field state is \xp). Because of 
Eq. 24, it follows that 

p( a | \i >) ) > 0, for a € C, (26) 


and 



ap(a ||^)) = 1, 


(27) 


Here dj and Ij , for j = 5, / , are the annihi- 
lation and identity operators for the signal 
mode (frequency y), and the image mode 
(frequency v — 2i'//r)> respectively. Both 
of these modes beat with the LO to pro- 
duce IF waveforms, just as is the case in 
classical superheterodyne radio reception. 
It is easily verified that the real and imagi- 
nary parts of y are commuting observables 
on the joint state space, Hs ® 'Hr, and 
so are simultaneously measurable in the 
usual sense. Ordinarily, only the signal 
mode carries information, i.e., the image 
mode is unexcited. Under these circum- 
stances, the PDF for the observed y- value 
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reduces to 8, 9 

pCi/IIVOs) = -IsfalVOsI 2 . 

7 T 

for y € C, (29) 

where |y)s is the signal- mode coherent 
state with eigenvalue y, and |t/>)s is an 
arbitrary signal-mode state. Comparison 
of Eqs. 29 and 25 completes the demon- 
stration that heterodyning — a pair of com- 
muting observables on an extended state 
space with an unexcited image mode — 
realizes the ds-POM. 

3* PHASE ESTIMATION 


can be taken to be { dH(<f>) : — 7r < <f> < 
7 r } , where 

= dn(^») t , (32) 

and 

/ = r rfnw, (33) 

J — 7T 

on the state space of the output mode, a. 
The conditional probability density, given 
$, for obtaining a phase value <f> from this 
POM is 

for -it < <£, $ < 7r, (34) 


Rather than seeking a quantum for- 
malism for phase by pursuing a quantum 
version of Ekp 1, Shapiro, Shepard, and 
Wong 3, 4 approached the problem from 
the estimation theory viewpoint. Con- 
sider the following abstract quantum esti- 
mation problem. A single-mode input field 
of annihilation operator ajN and quantum 
state undergoes an unknown, non- 

random, c-number phase shift $>, yielding 
a single-mode output field of annihilation 
operator 

a ~ e l ^d//v, (30) 

in state 

|V>) = exp (i$a\ N aiN} \ ip) in- (31) 

By making an appropriate quantum mea- 
surement on the a mode, and knowing the 
input state (V 7 ) 7 at, we are to estimate the 
phase shift <F. An interferometric phase 
measurement can be embedded into this 
scheme by placing appropriate constraints 
on the allowable quantum measurement. 
Optimizing a phase measurement within 
this more restricted environment cannot 
outperform the behavior obtained from 
an unfettered measurement optimization. 
Indeed, we should expect that joint op- 
timization of the quantum measurement 
and the input state will yield superior 
phase estimation performance. 

Without loss of generality, we can con- 
fine the phase shift to a 27r-rad interval, 
i.e., we can assume that — it < $ < it. The 
class of POM’s we must optimize over, in 
order to find the best phase measurement, 


where |^) is the state of the a-mode. 

In classical estimation theory, the 
maximum-likelihood (ML) estimate $ ml 
of an unknown, nonrandom, phase shift <!>, 
based on a noisy phase-shift observation 0, 
of known PDF p( 0 | $ ), is the phase shift 
which maximizes the likelihood of getting 
the observed datum, i.e., 

= arg max p(<j> \0). (35) 

— n<0<TT 


Often, the ML phase estimate equals the 
observed phase shift, because p(</>\ 9) has 
its peak at 9 = <f> } for —it < <j> < it. Such 
is the case for phase estimation in additive 
white Gaussian noise. 11 It then follows the 
p(4> | <t>) } the peak likelihood, is a simple, 
but meaningful, performance measure for 
$ M L* Indeed, its reciprocal, 

- p{<t>\ <f>) 

is the PDF’s width for the case of a uni- 
form distribution; if the distribution is 
Gaussian, then we have 6<fr = \/27rA </>, 
where A0 is the root-mean-square (RMS) 
error. 

Our problem is one of quantum estima- 
tion theory, namely, choosing the POM, 

<ill(0), and the input state, \rp)iN, to opti- 
mize our estimate of the phase shift For 
a given POM and input state, Eq. 34 sup- 
plies the PDF needed to perform classical 
ML estimation. In this quantum setting, 
however, the observed phase value 0 is, by 
presumption, our estimate of <I>. Thus, in 
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order for this estimate to be one of maxi- 
mum likelihood, we can restrict our atten- 
tion to POM’s satisfying 

<j> = arg max p(<j> \ 0), 

— w<0<ir 

for — 7T < <t> < 7T, (37) 

and optimize our estimate over dll 
and | ip) in by maximizing the peak 
likelihood — minimizing 6<p — averaged over 
all possible $ values. Here it is known 
that, for the input state whose number 
representation, xp n = (n|t/>), is 

rp n = |V>n|c ix ", fom = 0,1,2, , (38) 

6(f> is minimized by the following POM/ 

dn(<t>) = le^Xe*,^, 

for — 7T < 4> < 7T, (39) 

where 

|e^,^> = f> <(n * +XB) |n>. (40) 

n=0 

Moreover, the reciprocal peak-likelihood 
that results when we use this optimum 
POM to estimate $ is easily shown to be 

6<j> = 27r|(e"*’,i/#)| 



which is independent of the phases {xn}- 
In fact, p(<p | $) is independent of the 

{Xn}- 

We can exploit the {xn} independence 
to good purposes by assuming, without 
loss of generality, that the input state has 
positive real ip n . Equation 40 then reduces 
to 

OO 

I = |e‘») 2 •£ 

n—0 

for — 7r < <t> < ir, (42) 

which is the number-ket expansion of 
the SG phase operator’s (infinite-energy) 
eigenkets, viz. 

e^|e^) = 

for —it < <f> < it. (43) 


This says that the SG-POM is the quan- 
tum measurement for ML phase estima- 
tion in the general measurement configu- 
ration when the input state has a posi- 
tive real number representation. In other 
words, the phase eigenkets of the SG oper- 
ator generate the resolution of the identity, 

f = 7 L f d^Xe^l, (44) 

needed for ML quantum phase estimation 
in this case. For an arbitrary input state, 
the optimum POM from Eq. 39 is equiva- 
lent to performing the unitary state trans- 
formation 

OO 

U = e - * x "|n)(n|, (45) 

n=0 

followed by the SG-POM. 

To achieve the goal of jointly optimizing 
phase-estimation performance over both 
the measurement and the input state, it 
only remains for us to minimize 64>, from 
Eq. 41, by appropriate choice of \ip)iN- 
This problem has been addressed/’ 4 and 
the state 

A 


for n = 0, 1 , 2 M < oo, (46) 

where A is a normalization constant and 
M is a truncation parameter, has been 
shown to achieve 

6<p ~ 1 / dV 2 , (47) 

in terms of its average photon number, 
N — (a* a). This performance is far supe- 
rior to the 6<f> ~ \/N reciprocal peak like- 
lihood capability of optimized squeezed- 
state interferometry. However, the phase 
measurement PDF for the Eq. 46 state 
is a heavy-tailed distribution, viz., its 
RMS phase error, A<£, is essentially in- 
dependent of N. Thus, the degree to 
which this reciprocal peak likelihood ad- 
vantage can be usefully exploited has 
yet to be established. 12 In what follows, 
therefore, we will concentrate on the SG- 
POM, in that it constitutes the maximum- 
likelihood quantum phase measurement 
for all quantum states. 


Ill 


4. SG POM 


resolve the identity is also easily shown, 


The Susskind-Glogower operator, 2, 13 

e*^ = (od^) -1//2 d, (48) 

affords a well-defined polar decomposition 
of a, 

d = + (49) 

in terms of energy (number) and phase op- 
erators. Using the number-ket expansions 
of N and a we have that 


r d* |e*)(e*| 

J —7V 

oo oo 

= EE 

n=0ra=0 J '’ 

OO 

= 27r ^2 |n)(n| = 2*7. (57) 

n— 0 

That they are not orthonormal can be 
demonstrated from some simple Fourier 
transform manipulations, 14 


|n)(n + l|, 

n=0 

from which it follows that 
e~ x * = e'* 


(50) 


= £ e 
n— 0 

1 




= X! I n + ^ ei<t> > ( 51 ) 

n— 0 


= 

\n — — oo 

oo > 

+ sgn(n)e _,Tl ^ + 1 

n= — oo / 

= *■£(<// — </») 


and 




jot ('*':*) 

= |0)(0|. 

(52) 

2 V 2 ) 


(58) 


In words, the SG phase operator is not Here, ^(‘) is the Dirac delta function, and 
Hermitian, and does not commute with its 

adjoint. Thus, as was seen earlier for d it- f — 7 n < 0, 


self, the quadrature components of the SG sgn(n) 

operator, 


f -1, n < 0, 
= < 0, n = 0, 
l 1, n > 0, 


(59) 


cos(4>) = Re(e*^), and (53) 
sin(<£) = Im(e*^), (54) 

are noncommuting observables, 

[cos(^),sin(</>)] = ||0)(0|. (55) 

^The SG-POM derives from the fact that 
e 1 ^ has an overcomplete set of eigenkets, 
cf. the Sect. 2 discussion of the a-POM. 
By direct substitution of Eq . 50, we can 
verify that 

e^le*) = e^|e^), 

for — 7 r < <t> < 7r, (56) 


is the signum function. 


5. COMMUTING OBSERVABLES 

Recall that a POM on H which is not 
an observable on that space can be repre- 
sented as a collection of commuting ob- 
servables on a larger, signal x apparatus 
state space, H ® Ha> with the apparatus 
placed in some appropriate state. We now 
develop such a representation for the SG- 
POM. Aside from alleviating the qualms 
of those who believe only in observables, 
this representation may guide us to a real- 
ization of the SG-POM — the commuting- 
observables description of the a-POM is 
intimately connected with its heterodyne- 
detection realization. 


where |e l<i,!> ) has the number-ket represen- Let oa be the annihilation operator of 
tation given in Eq. 42. That these kets an apparatus mode, whose state space, 
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Ha, is spanned by its number kets, { |n)^ : where V'nm = /t(m|(n|i^5 X yv)- These two 

n = 0,1,2,...}. The non-Hermitian oper- are properly normalized in that tor all 
ator | ip) s xA we have that 


Y = (e*^ ® |0) i 4A{0|) 

+ (|0)(0| ® e~^A^j , (60) 

where 

= a\{aAa\)~ x/2 , (61) 

is easily shown to commute with its ad- 
joint. Here, Y is an operator on the joint 

state space H ® "Ha, and e'^A is the ap- 
paratus mode’s SG phase operator. 

Because [y,yt] = 0, the quadrature 

components of Y — denoted Y\ and Y^ 
are commuting observables, which can be 

measured simultaneously, i.e. , Y — Y\ +Y% 
can be measured in the usual sense. Solv- 
ing for the eigenkets and eigenvalues of Y 
we find that signal x apparatus number ket 

|y) = |n)|m)/i, for nm > 0, (62) 

is a T-eigenket with zero eigenvalue, and 

1 


\Y) = 




(|0)|0)A 


Pr(0 | Wsxa) 

+ f d<l>p(<f> | \i>)sxA ) = 1, (66) 

J — 7T 

as required by classical probability theory. 

Now, the commuting-observables repre- 
sentation of the SG-POM is at hand. Sup- 
pose we measure Y when the apparatus 
mode is unexcited, i.e., \ip)sxA = |V , )|0)>ii 
where | %p) € H is an arbitrary signal-mode 
state and |0),t is the apparatus mode’s 
vacuum state. Then the discrete value 
zero is never obtained, and the PDF for 
obtaining Y = e l<t> reduces to 

p(<j) | |V>)|0)a) 

= s|£ «-*■**• 

= 

for —7 r < 4> < 7T, (67) 


+ E ( ein ^l«)l°)A + e-^|0)|n) A ) ) , 

n=l ' 

for — 7r < 4> < 7T) (63) 

is a T-eigenket whose associated eigen- 
value is . Collectively, these com- 
prise a CON set from which we have 

that measurement of Y, when the 
signal x apparatus state is | il>)sxA G H ® 
Ha, yields a mixed classical random vari- 
able, y, which takes on either the discrete 
value 0, or a value from the continuum 
{ e*^ : — 7T < <p < ir }. The former occurs 
with discrete probability 

OO OO 

Pr( 0 I WsxA )=EE l^l 2 : C 64 ) 

n— 1 m=l 

the probability density for the latter is 
p(<t> I Ms* a) = ^ lV»oo+ 

E (e- in ^n0 + e <n ^0n) 


n = l 

for —7 r < <p < 7r, 


realizing the SG-POM statistics for an ar- 
bitrary state of the d-mode. 

The equivalence of the SG-POM to the 
Y measurement with an unexcited appar- 
atus mode allows us to clarify some ba- 
sic points. First, because the SG oper- 
ator does not commute with its adjoint, 
it is really the operator analog of the c- 
number e 1 ^ from the classical single-mode 
field. In other words, there is no Hermi- 

tian phase operator, 0, on such that 
exp(i<£) = (dd t ) _1/2 a. Restated in terms 
of the quadratures of the SG operato^this 

means that cos(0) 7^ cos(0), and sin (0) 7^ 
sin (<j>). As a result, the classical trigono- 
metric identity, 

cos(<£) 2 + sin (4>) 2 = 1, 

for — 7T < (f> < 7T, (68) 

does not apply to the quadratures of the 
SG phase operator, e.g., because 


(65) 


-7^2 ~ 2 |0)(0| 

cos (<f>) + sin (<£) = I — 


(69) 
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any \rp) with 7 ^ 0 gives 

(^|cos(0) I V>) + (}p\sm(<p) 2 \rp) < 1 . (70) 

However, the outcome of the SG-POM is 
a phasor e 1 ^* — the Y measurement with 
a vacuum-state apparatus mode yields a 
complex- valued, continuous classical ran- 
dom variable V = e 1 ^, where — 7 r < <f> < tt. 
Thus, we have that \Y\ 2 = 1 , with proba- 
bility one. 

The second point to note regarding the 
SG-POM is the relation of its mean value 
to those of the SG quadratures. Using 

U = (cos(<£) <8 |0 )m/i(0|) 

+ (| 0 )( 0 | ® cos( 0 )^ , (71) 

and 

Y 2 = (sin(0) <8> |0) , 4 . 4 ( 01 ) 

- (|0)<0| (S> sin(0)^) , (72) 

and assuming an unexcited apparatus 
mode, we find that 

/l(0|(V'|i > i|V’)|0)/i = <-0|cos(^){-0>, (73) 

and 

A(0MY 2 m0)A = (i/'|sin(^>)|V>), (74) 

for all o-mode states, 1*0) . What this says 
is that averages of the classical cos ((f)) and 
sin(0) random variables obtained from the 
SG-POM coincide with averages of the 
SG phase operators quadratures. To the 
extent that the quadrature mean values 
comprise the information of interest, we 
can conclude that the SG-POM provides a 
proper quantum measurement description 
for simultaneous extraction of this infor- 
mation from both quadratures. 


6. UNCERTAINTY PRINCIPLE 

Number-ket expansions of the quadra- 
ture operators cos(^) and sin(0) lead to 
the commutator 

[n , sin(<£)J = ico6(0), (75) 


and the associated uncertainty principle 

(AN 2 )(Asm(<p) ) > i |(cos(<£))| . (76) 

Equation 76 is valid for arbitrary states, 
but its utility, in this general form, is 
somewhat limited. First, the minimum 
uncertainty product is state dependent — 
a consequence of Eq. 75 not being a o 
number commutator. Second, the princi- 
ple does not directly address the variance 

of a phase measurement — it is the sin(0) 
operator whcee variance appears. 

It is common practice to use the lin- 
earized form of Eq. 76, 

ANA<t> > i , (77) 

which applies for states meeting the high- 
mean-field condition, 

(N) « |(a)| 2 » 1. (78) 

The linearized result, while useful, can be 
abused. Number kets have zero number- 
measurement uncertainty, and (a) = 0, 

(cos(<£)) = 0, hence the general result 
leads to the correct number-ket limit, 

(AlV 2 )(Asin(^>) 2 ) > 0, (79) 

whereas the linearized form is inapplica^ 
ble. 

Although the SG-POM does not allevi- 
ate the state-dependent nature of number- 
phase uncertainty limits, it does lead to 
an uncertainty principle which directly 
addresses phase variance. Our route to 
this principle — through Fourier theory — 
has the following motivation. The 
time-bandwidth uncertainty principle for 
the continuous- time Fourier transform 
(CTFT) 14 can be applied to the normal- 
ized position and momentum wave func- 
tions, i^(ai), and ^( 02 ), because they sat- 
isfy the Fourier transform relations 5 

*(a 2 ) = 

f°° ^rP(cn)e-' 2a ' a >, (80) 

J —00 v* 

and 

V»(ai) = 

/ ^*(a 2 )e i2aiQ2 . (81) 

J — 00 V 7T 
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The result of this procedure can be re- 
duced to 


(Ao?)(Ad^) > 


(82) 


= S>-n) 2 hM 2 

=0 

: f ^(*-*) 2 |*(e*)| a ( 87 ) 

J-i f27T 


n=0 
X 


which is the Heisenberg uncertainty 
principle for the annihilation operators 
quadratures. 

Because of Eq. 44, any state \ijj) has a 
phase representation 

^(e^) = 

for — 7T < <p — ^1 

such that 

IV») = W) 

= — r dtn**) i e< *>- t 84 ) 

27T J-it 

The phase representation of \ip) is inti- 
mately related to its number-ket represen- 
tation, = (n|t/>) they are a Fourier 
transform pair 


r d$ IdVie'*) 

J- w 2 tt 


d4> 


x ( 8g ) 

J- % 2n 


d4> 


(89) 


> ^]}’ 


2tt 


Ip( tt I |i^> > — I] 2 . 


(91) 


— irup 


(85) 


tf(e^) = 51 ^" e " 

n=0 

and 

i /" <#*(e‘V* (86) 

Z7T 7 — 7r 

as can be seen from Eqs. 83 and 42. In 
other words, ^(e 1 ^) and V'n constitute 
phase and number wave Junctions, which 
are capable of representing arbitrary 
states. The complementarity of the num- 
ber operator measurement — whose proba- 
bility distribution is Pr (N = n) = |V>n| 2 
and the SG POM— whose probability den- 
sity function is p(</>) = |'l'(e“ # ’)| 2 /27r— then 
follows from the Fourier relations, Eq. 85 
and 86. Thus, to obtain a number-phase 
uncertainty principle for the product of 
the number-operator variance and the SG- 
POM variance, we shall exploit this com- 
plementarity by paralleling the standard 
Fourier derivation of Eq. 82. 

With (AN 2 ) denoting the number- 
measurement variance and (A 4> 2 ) the SG- 
POM variance, when the field is in an ar- 
bitrary state | i>), we have that 

(AN 2 )(A<t> 2 ) 


In this development: n and <f > 'are the 

mean values of the number and SG-POM 
measurements, respectively, on the state 
uy ^'(e^) = the Schwarz 

inequality has been used in Eq. 89; and 
the integration necessary to obtain Eq. 91 
follows from the SG-POM’s PDF, p( 4> | 
|t/,)) = |*(e^)| 2 /2*- 

Unlike the usual number-phase uncer- 
tainty principle, i.e., Eq. 76, our result 
does not require any linearization before 
it can be applied to phase variance. Equa- 
tion 91 is still state dependent, but this is 
unavoidable. When \rp) is a number state, 
we have 

p( <t> I M ) = 2^> for ” 7r < ^ - 7r ’ ( 92 ) 

a uniform distribution, which is maxi- 
mally random, but still has finite variance. 
There is no contradiction with Eq. 91 in 
this case, even though AN = 0 for a num- 
ber ket; the uniform PDF causes the right 
member of Eq. 91 to vanish. On the other 
hand, when |t/>) is a high-mean-field state, 


we will have p(ir 
Eq. 91 reproduces t 
formula, Eq. 77. 


. |V>>) « 1, so that 
he standard linearized 
Indeed, for any state 


iormuia, dm. — r * . . 

satisfying p( JT |J) « 1, we have that 
Eq. 77 holds. This makes the SG-POM 
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derivation of Eq. 77 more robust than lin- the connection with causal waveforms, we 
earization of Eq. 76 . introduce the term number-ket causality 

for the condition Eq. 94. This is a well- 
studied problem in linear systems, so re- 
7. Number-Ket Causality sults are immediately available. 14 ’ 15 


The SG-POM underlies maximum- 
likelihood quantum phase measurement 
for all quantum states. Given the prob- 
lems associated with minimizing the SG- 
POM’s reciprocal peak-likelihood 6<fi, by 
choice of input state, 3, 4> 12 a different 
state-selection criterion may be worth con- 
sidering. In this vein, it is germane to ask 
the following question. What SG-POM 
phase PDF’s can be realized by choice of 
input state \ip)? It turns out that linear 
system theory has the answer. 


From Eq. 93 and the Paley- Wiener 
theorem 15 we have that p( <p | |^) ) must 
satisfy 

/ I ln[p( 4 >\ IV') )]| < oo, (95) 

for all number-ket causal ^(e^). From 
this condition, it follows that no state can 
confine the phase-measurement PDF to a 
subinterval of (— 7r, 7 t], e.g., the uniform 
density, 


The Susskind-Glogower probability op- 
erator measurement on a state |t p) results 
in a classical random variable with prob- 
ability density function 


p { 4 >\ h/>)) = 

for \4>\ < 6<j>/ 2 < 7 r, 


(96) 


p ( 4 > IM) = 


|*Ceg)g 

2n 


for -it < <f> < 7r. (93) 

Here, is the phase representation 

of the state j^). According to Eq. 86, the 
phase representation is the Fourier trans- 
form of the number representation. The 
latter is a one-sided , discrete-parameter 
sequence that is the inverse Fourier trans- 
form of i.e., 


is impossible . The Paley-Wiener condition 
is both necessary and sufficient, i.e., if a 
PDF obeys Eq. 95, then there is a state 
which gives this density through Eq. 93. 
Indeed, there are an infinite number of 
such states, because Eq. 93 only constrains 
the magnitude of the phase representation. 
One such state can be obtained explicitly 
via the discrete-parameter Hilbert trans- 
form. The procedure is as follows. For a 
PDF obeying Eq. 95, set 


V'n = Ty" 

(n|V>), for integer n > 0, 

(94) 

0, for integer n < 0. 

In system-theory parlance, the Fourier 
pair V’n 4 — ♦ ^(e 1 ^) is analogous to 
that for a discrete-time waveform on 
an unbounded interval and its periodic, 
continuous-frequency Fourier transform, 
i.e., the discrete-time Fourier transform 
(DTFT). 14 More importantly, saying that 
r/j n is one-sided is equivalent to saying that 
a discrete-time waveform is causal, viz. it 
could be the impulse response of a causal, 
linear time-invariant system. Determin- 
ing what p( <p | |^>) ) are possible from 
Eq. 93 is then the same as determining 
what |'J>(e'^)| are Fourier-transform mag- 
nitudes of one-sided {V’n}- To emphasize 


for -n < <j> < 7r. (97) 

Next, find the discrete-parameter Hilbert 
transform, 14 

arg[^ , (e*^)] = V j* ^ ln[|*(e*')|] 

xco t (^i), 

for -7r < 4> < ir, (98) 

where V denotes Cauchy principal value. 
Equations 97 and 98 then comprise the 
magnitude and phase — the polar form — 
of a properly-normalized phase represen- 
tation 'k (e 1 ^) with the prescribed phase- 
measurement statistics and a number-ket 
causal inverse Fourier transform. 
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The preceding phase representation 
construction for a state with prescribed 
SG-POM statistics is by no means unique. 
Equation 97 constrains |'I'(e < ^)|, but no re- 
striction is placed on argf^e’^)]. Con- 
sider a number- ket causal function, { h n : 
n = 0, 1, 2, ... }, whose Fourier transform 
has unity magnitude, viz. 

OO 

H(e *+) s Yl h ” e ~ in *' 

n— 0 

for —it < (p < it, (99) 

obeys 


there is an uncountable infinity of states 
which have the same SG-POM statistics. 
Nevertheless, the state constructed via the 
discrete-parameter Hilbert transform has 
a unique advantage — it is the minimum 
average photon- number state with the pre- 
scribed phase-measurement PDF. 

The proof of the minimum average 
photon-number property follows almost 
immediately from available linear-system 
results. Let : —it < 4> < Tt } 

and {rp n : n = 0,1,2,...} be the phase 
and number representations of the state 
\xjj), obtained via Eqs. 97, 98, and 94, 
that realizes a particular phase PDF. Sim- 


|tf(e’*)| = 1, for — 7r < <p < it. (100) 

Such a function is known in digital-filter 
theory as an all-pass filter; were { h n : 
n = 0, 1, 2, ... I the impulse response of a 
discrete- time, linear, time- invariant filter, 
the associated frequency response would 
pass all frequencies with neither attenu- 
ation nor gain. The prototypical exam- 
ple of an all-pass filter is obtained — in the 
z-transform domain — by balancing H (z)- 
poles within the unit circle with H(z)- 
zeros outside the unit circle to achieve 14 

K e~ i{ * - d£ 

= nfi^ 

for -Tt < 4> < It, 


ilarly, let {*'(e*) : -it < <f> < it } and 
{ip' n : n = 0, 1, 2, ... } be the phase and 
number representations of any other state, 
\ip'), with the same SG-POM statistics. 

Then, we have that 14 


Af-l 

E (w 2 - Wl 2 ) > o, 

n=0 

for M = 1,2,3,... (103) 

Physically, this says that, of all states with 
the desired phase behavior, the Hilbert- 
transform generated state concentrates its 
number-ket content closest to the vacuum. 
Because both Up) and \xp') are normalized, 
: ~ Vn 103 is eauivar- 


where \pk\ < 1, 


lent to 


for fc = 1,2, . . . tf.(101) 

Now, suppose we assemble the phase rep- 
resentation 

*V*) = 'J '(e^)//(e^), 

for —it < <p < it, (102) 

where ^(e*^) is constructed according 
to Eqs. 97 and 98 for a desired phase- 
measurement PDF, and H{e *^) is an all- 
pass phase representation from Eq. 101. 
The convolution-multiplication theorem of 
Fourier analysis, plus the fact that con- 
volving two causal functions produces a 
causal function, 14 guarantees that ^'(e 1 *) 
is a properly normalized, number-ket 
causal phase representation; the all-pass 
nature of H(e i< ^) implies that 'P , (e' 1 ^) has 
the desired SG-POM statistics. Because 
this process holds for all K j> 1 and for 
all pole locations within the unit circle, 


A„ = E ( W 2 - W.I 2 ) < 0, 

n—M 

for M = 0,1,2,... (104) 

Proving the minimum average photon- 
number property is now straightforward: 

(4>\n\xp) - wm’) 

= E "(i^i 2 - W.I 2 ) 

n=0 

oo 

= £ Am < 0. (105) 

M— 1 

Thus, Eqs. 97 and 98 provide the means 
for choosing a state of minimum average 
energy and prescribed phase-measurement 

PDF. 
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8. PHASE COMMUNICATION 


Sections 1-7 constitute an abridged ver- 
sion of Shapiro and Shepard. 4 That par 
per presents additional details regard- 
ing the state that achieves 6(f> ~ 1 /TV 2 , 
as well as substantial material on new 
classes of quantum states — coherent phase 
states, squeezed phase states, rational 
phase states — that are closely associated 
with the SG-POM. Furthermore, it proves 
that the Pegg-Bamett Hermitian phase 
operator 16, 1 '—which exists on a trun- 
cated state space and provides phase- 
measurement statistics on the full state 
space through a limiting procedure — is 
included within the SG-POM formalism, 
i.e., these two schema produce identi- 
cal phase measurement statistics for all 
quantum states. Neither of these top- 
ics will be considered herein. Instead, 
we shall move away from the single- 
mode case and develop new results for 
two-mode quantum phase measurement. 

Our objective will be to exploit the Y - 
measurement — developed en route to the 
commuting observables form of the SG- 
POM — when the signal and apparatus 
modes are quantum-mechanically corre- 
lated and a phase-conjugate modulation is 
applied to them. 

Consider the phase-conjugate quantum 
measurement setup shown in Fig. 1. This 
is a phase-conjugate system because what- 
ever c-number phase shift $ is applied to 
the signal mode, leading to the annihila- 
tion operator transformation 


motivation is to circumvent the Paley- 
Wiener restriction that encumbers phase- 
measurement PDF’s for single-mode fields. 
As we shall see, some startling new possi- 
bilities arise with two modes. 

The Paley-Wiener condition applies to a 
single-mode phase PDF because this den- 
sity is proportional to the squared magni- 
tude of the Fourier transform, {^(e^) : 
—7r <<£<*•}, of a one-sided sequence, 
{tp n : n = 0, 1,2,..., }. We shall break 
out of this limit, in the two-mode case, 
through quantum correlation. On Hs ® 
Ha, the joint state space of the signal and 
apparatus input modes, we can construct 
number-product vacuum states of the form 

oo 

Win = V>o|0) 5 |0U + X> n |«>s|0>* 

n=l 

+ rp- n \0)s\n) A ) , (108) 

where 

oo 

£ N>nl 2 = i. ( 109 ) 

n=— oo 

The term number-product vacuum is ap- 
propriate for such | xJj)in because, when 
the signal x apparatus state is of this class, 
a measurement of the number-operator 
product — N$ ® Na — yields outcome zero 
with probability one. Thus, for \xp)iN 
a number-product vacuum state, Eqs. 62 

and 63 imply that measurement of Y 
yields a classical phasor e l< ^, with —n < 
<f> < 7T. Moreover, <f> in this case has PDF 


-IN 

a S 




(106) 


the conjugate phase shift, — $ is applied to 
the apparatus mode, viz. 


-IN 




(107) 


cf. Eq. 30. If we take the signal and 
apparatus modes to be the appropriate 
linear polarizations, a transverse electro- 
optic modulator can be used to induce 
the necessary conjugate phase shifts. 18 
Phase-conjugate shifts also appear, pro- 
totypically, in gravity-wave detecting in- 
terferometers. In fact, there are funda- 
mental advantages to operating a phase- 
sensing interferometer in phase-conjugate 
fashion. 19 Our work does not depend ex- 
plicitly on the means by which this mod- 
ulation is accomplished. Its principal 


pf«H4>) = ^ . 

for — IT < <f>, $ < 7T, (110) 


in terms of 

OO 

*(e‘*) = £ V> n e-’"*, 

n=— oo 

for -7T < (j) < 7T. (Ill) 

Note that {rj> n : \n\ = 0,1,2,...,} 
and : —n < <j) < tt} are 

not tne number and phase representa- 
tions, respectively, of any single-mode 
field state. They are, however, the number 
and phase representations, respectively, 
for a two-mode, number-product vacuum 



state. The {VVi> ^(e 1 ^)} notation is con- 
venient because, as shown by Eq. Ill and 
its inverse, 



for |n| = 0, 1, 2, ... , (112) 

these functions are a Fourier transform 
pair. More importantly, this notation 
makes clear the fact that number-ket 
causality does not restrict the possible 
two-mode phase PDF’s. In particular, 
there are number-product vacuum states 
that satisfy 



0, for \<p\ > 4 > c , 
with cf) c < 7r. (113) 


i.e., two-mode phase PDF’s can be con- 
fined to subintervals of (— 7r, i r], a situation 
that is forbidden to the single-mode case, 
cf. Eq. 96. This possibility is of great sig- 
nificance for phase-based digital communi- 
cation and phase- based precision measure- 
ment, as we shall see. 


To cast the Fig. 1 structure into a dig- 
ital communication mold, let us assume 
that transverse electro-optic modulation is 
used to transmit a randomly- selected digit 
k , satisfying 1 < k < if, by using 
where 


probability can be achieved at finite root- 
mean-square (RMS) photon number. 21 

In order to achieve zero error probability 
in phase-conjugate quantum communica- 
tion, we need only use a number-product 
vacuum state which enforces Eq. 96, with 
<t>c < n/K. Under this condition, we have 
that the observed phase, 0, satisfies 



We also have that \$k “ > 2 7r/if, for 

j ^ k . So, for any observed 0 in the inter- 
val (—7r.7r], we know that 

Pr($ = *j | < |r) = 0, 

for all j ^ k. (116) 

This means we can unambiguously deter- 
mine which digit was sent by choosing the 
index associated with the unique $-value 
that is within 7 r/ifrad of the observed 
phase. Via this procedure we decode k 
from <j> with zero probability of being in- 
correct. 

To make this technique for zero error 
probability communication more explicit, 
we shall introduce a specific input state 
which has the desired property. The num- 
ber representation we shall presume is 


= -7T + 


(2 k - 1)tt 
K 


(114) 


Our objective is to make a minimum er- 
ror probability decision as to which 
was sent, based on the result of the Y - 
measurement when the signal x apparatus 
state is a number-product vacuum, char- 
acterized by {rpn, ^(e 1 ^)}. 


Hall and Fuss have considered the 
single-mode version of this if- ary dig- 
ital communication problem 20 . They 
optimized a single-mode state — in con- 
junction with the SG-POM — to obtain 
a phase-based quantum communication 
setup whose error probability vs. average 
photon number is significantly better than 
that for optical heterodyne detection. Hall 
and Fuss found a nonzero error probability 
at finite average photon number, TV, which 
approached zero as iV — ► oo. Surpris- 
ingly, in our two-mode problem, zero error 



for |n| = 0, 1,2, , 
and if = 2,3,4,... (117) 

The associated phase representation for 
this state is easily computed to be 


*(e^) = 

f ^Fcos(^) , for \4.\ < % (iig) 

[ 0 , for 7 ^ < \ 4 >\ < 7r * 

In Fig. 2 we have plotted p( 0 | $ = 0 ) vs. 
<f> for this state when if = 4; we see that 
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the nonzero support of this PDF is the in- 
terval (— 7 t/ 4, 7 r/4). In Fig. 3 we indicate 
how this PDF leads to zero error probabil- 
ity phase-conjugate communication; this 
figure plots the four conditional PDF’s, 
{p(0|4>fc):l<A;<4}, which apply 
when K — 4. For any observed Rvalue 
we must have that p{4> | <J>) > 0, other- 
wise that <£- value could not have occurred. 
Figure 3 shows that, for any given <p y 
there is only one possible 4>fc-value which 
satisfies the nonzero PDF requirement — 
zero error probability communication re- 
sults from deciding that this value was the 
transmitted phase. 

The next question to address is the pho- 
ton number statistics associated with our 
phase-conjugate communication scheme. 
For the general number- product vacuum 
state we have that the total — signal plus 
apparatus — photon number measurement 
has the following probability distribution, 

Pr(AT 5 + N a = n ) = 

IV>o| 2 , for n = 0, 

(119) 

iV'-nl 2 + iV'nl 2 , for n = 1,2,3,... 

For the particular state given by Eq. 117, 
it is then a simple matter to show that 

{N s + N a ) < v^s + ^a) 2 ) 

= f- (120) 

In other words, we can achieve zero error 
probability K- ary phase-conjugate quan- 
tum communication with an RMS total 
photon number of K/2 . Figure 4 is a plot 
of Eq. 119 for the K = 4 case. 

The preceding quantum communica- 
tion result is, of course, idealized. We 
have presumed a state generator — to pro- 
duce a specific number- product vacuum 
state — that as yet has no explicit realiza- 
tion. Likewise, our scheme uses the Y- 
measurement; again, no explicit realiza- 
tion is yet available. At least we can say 
that electro-optic modulation will impress 
the phase information on the input state, 
once that state can be produced. On the 
other hand, we have implicitly assumed 
lossless transmission; inclusion of loss will 
inevitably lead to nonzero error probabil- 
ity. 


Our main purpose in going to the two- 
mode construct was to develop poten- 
tial quantum- phase measurement schemes 
that promise substantial benefits, i.e., ben- 
efits that warrant the effort to bring 
them to fruition. This motivation is 
very much in line with the starting point 
for squeezed-state research. 22 In this re- 
gard, it is instructive to compare our 
phase-based scheme for zero error prob- 
ability quantum communication with a 
more well-known approach based on num- 
ber kets. For a single-mode field with 
annihilation operator a, lossless transmis- 
sion of one of the number kets { \k — 1) : 
1 < k < K } followed by ideal direct de- 
tection, viz. the N = d^a measurement, 
also yields K- ary digital communication 
without error. For k equally likely to be 
any digit between 1 and the average- 
energy efficiency of such a single-mode, 
number-ket system is roughly the same 
as that of our two- mode, phase- conjugate 
system, i.e., both need slightly less than 
Kf 2 photons on average. The number- 
ket system has the advantage that its state 
generator may be approximated via feed- 
forward control using photon- twin beams, 
and its measurement only requires a high 
quantum-efficiency photon counter. Also, 
the number-ket approach uses less band- 
width; only one mode is needed. Alter- 
natively, number-ket direct detection on 
a two- mode field can be used for error- 
free K-&iy communication at significantly 
less than K/2 photons on average. How- 
ever, if we shift our attention from phase- 
based communication } to phase- based pre- 
cision measurements , the Fig. 1 arrange- 
ment has a capability that number kets 
cannot match — phase sensing with con- 
trolled precision. 

Suppose that we use the Fig. 1 ar- 
rangement for phase- conjugate precision 
measurement. Specifically, let us use the 
number- product vacuum state Eq. 117 in 
conjunction with a phase-conjugate in- 
terferometer (see, e.g., Bondurant and 

Shapiro 19 ) and the Y measurement. Now, 
the phase shift <f> takes on any value from 
the continuum (-7T, tt]. Nevertheless, ex- 
cept for 27r-modularity effects which come 
into play when <£ is within n/K of ±i r, 
the observed phase will lie within n/K rad 
of the true phase with probability one . 
Thus, using less than K/2 photons on av- 
erage, we can guarantee a phase measure- 



ment which is within ■n/ K rad of the exact 
value. In other words, unlike more con- 
ventional schemes — which only ensure an 
acceptable RMS phase-estimation error — 
our phase-conjugate interferometer pro- 
vides exact phase determination to a pre- 
scribed number of decimal places. 
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Figure 1: Phase-conjugate quantum communication system. 
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Figure 2: Conditional phase-measurement PDF, given $ = 0 , for the state Eq. 118 
when K = 4 . 
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Figure 3: Conditional phase-measurement PDF's for the state Eq . 118 when K - 4 
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Figure 4: Signal-plus-apparatus number distribution for the state Eq. 11 7 when K — 4 
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